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Abstract 

We present a formulation which allows heavy quark (c, b, . . .) mass effects to be ex- 
plicitly incorporated in both the coefficient functions and the splitting functions in the 
parton evolution equations. We obtain a consistent procedure for evolution through the 
threshold regions for cc and bb production in deep inelastic scattering, which allows the 
prediction of the charm and bottom quark densities. We use the new formulation to 
perform a next-to-leading order global parton analysis of deep inelastic and related hard 
scattering data. We find that the optimum fit has as(Af§) = 0.118. We give predictions 
for the charm components of the proton structure functions Fi and Fl as functions of 
x and Q 2 and, in particular, find that F| is in good agreement with the existing mea- 
surements. We examine the Q 2 range of validity of the photon- gluon fusion model for 
cc electroproduction. We emphasize the value of a precision measurement of the charm 
component F% at HERA. 



1. Introduction 

A very wide range of deep inelastic scattering structure function data can be successfully 
described in terms of universal quark and gluon distributions satisfying DGLAP (Q 2 ) evolution 
equations. While the formalism for light quarks (i.e. m q <C Aqcd) is on a sound theoretical 
footing, the treatment of heavy quarks (i.e. m q 3> Aqcd) is more problematic. The reason 
is that in practice one requires a consistent description which includes both the kinematical 
regions Q 2 ~ m 2 and Q 2 3> m 2 

The problem of how to treat heavy quark contributions to deep inelastic structure functions 
has been widely discussed, see for example M. It has been brought into focus recently by the 
very precise F 2 ep {x 1 Q 2 ) data from HERA. Both the HI and ZEUS collaborations have measured 
0, [| the charm quark component _F 2 C of the structure function at small x and have found it 
to be a large (approximately 25%) fraction of the total. This is in sharp contrast to what is 
found at large x, where typically F 2 C /F 2 ~ 0(1(T 2 ) [|. Since the HERA F 2 data g are 
a potentially valuable source of information on the gluon distribution, the value of as, and 
the relation between the non-perturbative (low Q 2 ) and perturbative (high Q 2 ) domains, it is 
important that charm component is treated correctly. 

In this paper we present a new, theoretically consistent method for calculating the heavy 
quark contributions to the deep inelastic electroproduction structure functions F 2 and F L .f\ 
Our main focus is on the charm quark contribution, although our results apply equally well 
for bottom and top quarks. The most important feature of our analysis is that it is applicable 
both to the threshold region Q 2 ~ 4m 2 , where phase space effects are important, and to the 
asymptotic region Q 2 3> 4m 2 , where the charm quark assumes the role of a massless parton 
and the DGLAP resummation of leading (a^lnQ 2 )" contributions is necessary. 

Before describing our formalism and presenting quantitative predictions, we briefly review 
existing techniques for treating the charm quark contribution to F 2 . The most simplistic ap- 
proach is to assume that a probe of virtuality Q 2 can resolve a charm quark pair in the proton 
sea when Q 2 > m 2 . Since such pairs originate from fluctuations of the gluon field, g — > cc, a 
perturbative treatment should be valid as long as m 2 3> Aq CD . As Q 2 increases, 0{m 2 /Q 2 ) 
corrections to the standard DGLAP evolution become less important, and the charm quark 
can be treated as a (fourth) massless quark. These ideas are embodied in the 'massless parton 
evolution' (MPE) approach 

c(x, Q 2 ) = for Q 2 </4, 

n f = 3 + #(Q 2 -/i 2 ) in P qg , P gg , ft, . . . , (1) 

where [i c = 0(m c ). The charm contribution to the structure function is then 

F^(x, Q 2 ) = lxc(x, Q 2 ) (2) 

1 Note that we only consider here the case of neutral current deep inelastic scattering. The case of charged 
current scattering, e.g. W~c — > s, will be somewhat different but can in principle be treated using the same 
techniques. 
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in lowest order. This is the approach adopted at NLO in the MRS (and CTEQ) global parton 
analyses, with /i c chosen to achieve a satisfactory description of the EMC _F 2 C data W. For 
example, in the MRS (A) analysis |7| it was found that /i 2 = 2.7 GeV 2 and that this was to a 
good approximation equivalent to taking 

2c(x,Q 2 ) = 6S(x,Ql) (3) 

with 5 ~ 0.02 at the input scale Q% = 4 GeV 2 . That is at the input scale, charm (c + c) 
was found to have approximately the same shape as the total quark sea distribution S, and 
moreover to form about 2% of its magnitude. The input parameter /i 2 (or equivalently 5) was 
chosen to give a good description of the EMC F% data. 

Although phenomenologically successful, the MPE model clearly cannot give a precise de- 
scription of the charm contribution in the threshold region. Two-body kinematics imply that 
an on-shell cc pair can be created by photon-gluon fusion (PGF) provided 

M /2 = g 2 l^> 4m 2 (4) 

where W is the 7* g — > cc centre-of-mass energy. That is, at small x, cc production is not 
forbidden even at small Q 2 < /x 2 where the MPE approach gives zero|. In the PGF approach, 
which was used, for example, in refs. |j| PH, F$ is calculated using the exact matrix elements 



and phase space for the process 7*g — > cc. In leading order in as we have 

F 2 c (x, Q 2 ) = J dz C g (z, Q 2 , fi 2 ) X - g g, ^ . (5) 

Note that the scale /i 2 at which the gluon distribution and the coupling as (in the coefficient 
function C g ) are evaluated is not specified at leading order, but one might guess that fi 2 = 0(m 2 ) 
is appropriate. We discuss a reasonable choice of /1 2 in more detail in section 3 together with 
the effects coming from next-to-leading order (NLO) corrections. In contrast to the situation 
for massless quarks, there is no collinear divergence in the leading-order 7*^ — > cc calculation: 
the integral over the transverse momentum of the produced cc pair is regulated by the quark 
mass: / dk^k 2 -, / '{k^ + rn 2 ) 2 . However, this in turn means that at very high Q 2 the leading- 
order contribution behaves as F% ~ as{^ 2 ) \n(Q 2 /m 2 ). Higher-order corrections also behave 
as (a$ ln(Q 2 / m 2 )) n , and fixed-order perturbation theory breaks down. In fact these large 
logarithms are precisely those which are resummed by the DGLAP evolution equations. Thus 
at large Q 2 we have to include the charm quark as a parton in DGLAP evolution. The exact 
next-to-leading order corrections to the PGF structure function are known |12|. Indeed very 
recently |R| this leading-twist analysis has been used to perform (as \n(Q 2 / m 2 )) n resummation 



for Q 2 ^> m 2 . However such an approach is of course not applicable in the threshold region 
Q 2 > m 2 , which we also wish to study because of the HERA data for F2 in this domain. 

2 In ref. || the MPE model was modified by the introduction of a smooth 'smearing' function which gave a 
gradual onset of the charm distribution from a low input scale Qq, namely Qq = 1 GeV 2 . 
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Our goal is to include the charm quark in parton evolution in a consistent way. First, in sec- 
tion 2, we discuss how to include the heavy quark mass in the Altarelli-Parisi splitting kernels 
in such a way as not to destroy the original parton interpretation, that is to ensure energy- 
momentum conservation etc. In order to be consistent with the MS factorization scheme the 
threshold for the onset of the charm distribution is Q 2 = m 2 . Although in fact we will see that 
charm partons cannot be resolved in the timescale (At ~ 1/AE ~ 2E/Q 2 ) of the fluctuations 
of the gluon into a cc pair until Q 2 > Am 2 . In section 3 we discuss the coefficient functions. 
The PGF contribution will be included in the coefficient function for the gluon distribution. 
Thus below the resolution threshold of the charm distribution, Q 2 < Am 2 , our result for cc 
production will not be zero but will agree with the PGF approach. However, at large Q 2 , as 



was noted in ref. |p^| , part of the PGF cross section is automatically generated by the evolu- 
tion of the charm distribution. To avoid double counting we must therefore subtract from the 
coefficient function given by PGF the contribution which is generated by evolution in this way. 
As a consequence, above the charm threshold a smaller and smaller fraction of _F 2 C will come 
from the direct photon-gluon fusion mechanism, and instead the main part will be generated 
by conventional parton evolution. In section 4 we use the new formulation to perform a NLO 
global analysis of deep inelastic and related hard scattering data. We find an excellent overall 
description with, in particular, a significant charm component of F2 in the HERA regime. The 
analysis allows us to predict universal charm and bottom quark distributions, c(x, Q 2 ) and 
b(x, Q 2 ). In section 5 we present the partonic decomposition of F£ as a function of Q 2 and, for 
completeness, compare the PGF model estimates. We also give predictions for F\. In section 
6 we study the charm component of the longitudinal structure function F^. Finally, in section 
7, we give our conclusions. 



2. The effects of the charm mass on evolution 

As mentioned above, our aim is to develop the appropriate formalism to describe deep inelas- 
tic scattering which incorporates the production of a heavy quark pair (which for definiteness 
we take to be cc) and which allows a universal charm parton distribution to be obtained from 
an analysis of these and other data. We can identify the charm mass effects in the structure 
functions F% L which describe such scattering from the following subset of integrations^ 



rik 2 r rlh 2 h 2 r rlk 2 



h 2 / (h 2 -I- m 2 \ 2 I h 2 



(6) 



Ti+1 



where kxi are the transverse momenta of the t channel partons. The mass of the charm quark 
enters in the k\ { integration which results from the g — > cc transition, see Fig. 1. For the 
example of the parton chain shown in Fig. 1 it appears that m 2 should also have been retained 
in the integration over k\ i+v However, we show below that this is only needed at next-to-next- 
to- leading order (NNLO) in a$- 

3 The general structure of the integrands is fc^/(propagator) 2 , which for massless partons ~ \jk\. The fcy 
in the numerator arises from the spin structure of the parton vertex. 
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First we recall the kinematic regime responsible for the leading-order (LO) result. LO evo- 
lution corresponds to the resummation of the leading logarithm terms, (aslnQ 2 ) n , which arise 
when the n emitted partons have strongly ordered transverse momenta (. . . -C k^^i <C k^ <C 
kj<i+i ■ ■ ■)■ If two of the partons were to have comparable transverse momenta, krj ~ kxj+x, 
then we would lose a \nQ 2 and obtain instead a NLO contribution of the form asias^Q 2 ) 111 - 
We may write the Altarelli-Parisi splitting functions 

Pa 'T • ",/>;,' : + ... (?) 

where P^ is the LO form and P^ gives the NLO correction. P^ includes virtual corrections 
to the vertex and propagators as well as the possibility of producing a second l s channel' parton 
with comparable transverse momentum. 

2.1. LO evolution with charm 

On the scale of the Altarelli-Parisi evolution in InQ 2 , we see that to a good approximation 

that is the presence of the charm mass simply cuts out the contribution from the region 
k\i J$ rn 2 c . Indeed we will find, at LO accuracy, that we have massless three flavour evolution 
for Q 2 < m 2 and massless four flavour evolution for Q 2 > m 2 c \ that is due to strong ordering 
(&Ti+i ^ ^Ti) we can neglect the charm mass in the k\ i+1 integration of Fig. 1. Therefore at 
LO the singlet evolution equations have the symbolic form 

g = Pgg® 9 + ^ Pgq® q + Pgc® C 

q 

q = P qg ®g + P qq ®q (9) 

C = Peg ® 9 + Pec ® C 

where q = u,d,s denotes the light quark density functions and c the charm density. We have 
abbreviated P<°) by P and / = (2-k / a s )df / dlnQ 2 . At LO the quark mass effects are simply 
encapsulated by 

P ci = P ci (m c = 0) 6(Q 2 - m 2 c ) (10) 
with % = g or c, and similarly for P gc . Also the virtual contribution to P gg must be modified 

Pgg = ••• -\n f S{\-z) (11) 

with rif = 3 + 6(Q 2 — m 2 ), and, of course, we must allow for the increase in the number of 
active flavours nj in the running of as- 

Although we show in ([|) only the equation for c, we note that each heavy quark (c, b, . . .) 



requires a separate singlet evolution equation [15] since their splitting functions have different 
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6 function contributions. 



2.2. NLO evolution incorporating the charm mass 

At NLO the inclusion of quark mass effects is a bit more complicated, although it turns out 
that we only have to take m c into account in P cg and then only in the LO part (Of course 
as a consequence we must adjust the virtual corrections to P gg )- The argument is as follows. 

We have to improve on approximation (§) of the k\ { integration in (|). To do this we divide 
the integral into two parts 

dklik 2 Ti r d{k 2 Ti + m 2 c ) r m 2 c dk 2 Ti 



(k 2 Tl + m 2 ) 2 J (k 2 Tl + m 2 c ) J (k 2 Tl + m 2 ) 2 



(12) 



where the first term gives the leading logarithm contribution that we discussed in section 2.1. 
To be specific we have 

2 d(k 2 Tl + m 2 c ) , Q 2 



for fcf^i <C m 2 c1 which is equivalent to the threshold factor 9(k^ i — m 2 ,) of (H). The second term 



in fll2l) , which is concentrated in the region k\ { ~ m 2 , gives a constant contribution. That is, 
it is a NLO contribution (containing a factor as without an accompanying InQ 2 ). It means 
that the m 2 c effects need only be evaluated in the LO (one-loop) part of the g — > cc splitting 
function, P^ . For instance consider the integration over fc-ri+i of Fig. 1 and the possibility of 
m 2 effects in P cc . Clearly if k^ i+l 3> m 2 c then the mass terms vr? c jk\ i+1 can be neglected. If, 
on the other hand, k^ i+1 ~ m 2 c then either k\ { -C m 2 c and c(x, k\^) = or k\ i ~ m 2 c and we 
lose two InQ 2 factors so that the contribution is NNLO, which we omit here. That is, at NLO 
there are no m c effects in P cc . A similar argument shows that this is also true for P gc . 

In summary, we have shown that at NLO P cc and P gc remain as in section 2.1, whereas 

P cg = P c W(m c ) + a s P£(m c = 0) 6{Q 2 - m\). (14) 

That is we need only evaluate the effect of the charm mass on the LO part of P cg . As a 
consequence of the change in P cg , we have to adjust the virtual correction to P gg by an amount 
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-8(1 -z) £ dzz(p£\z,m c ) - pW(«,m B =0)), (15) 



see section 2.3. This adjustment also restores energy-momentum conservation. 

We note that instead of a charm density based on the Renormalization Group (RG) equa- 
tions and leading-twist contributions^ we have introduced an arguably more physical parton 
density based on the leading ln Q 2 decomposition of the Feynman diagrams retaining full mass 

4 The charm density in the conventional RG approach has been discussed recently in ref. p3| , where the 
(leading twist) coefficients have been fully calculated at NNLO in the limit Q 2 3> mj,. 
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effects. This charm density is universal^ and the partonic momentum sum rule is satisfied. Our 
definition of the parton density coincides^] with the conventional (massless) RG-based definition 
for Q 2 3> m 2 c . 

We have seen that in our approach at NLO we need only consider m c ^ effects in the 
LO diagrams. It is straightforward to extend the formalism to allow for charm mass effects in 
NNLO evolution. We need only consider m c ^ effects in the NLO diagrams. That is we need 
to evaluate the "blocks" gg — > gg, gg — > cc, qq — > cc to 0{a 2 s ) with m 2 c included explicitly, but 
only in the region k\ i ~ m 2 . For example, for gg — > gg we would need to evaluate the diagrams 
shown in Fig. 2. 



2.3. Evaluation of quark mass effects in P c 



We note that heavy quark mass effects were studied in refs. [To, 1C] in terms of the anomalous 
dimensions of the moments of structure functions. However, it is difficult to apply the results 
to parton evolution, since in these early studies the mass correction plays the role of a higher 
twist contribution. As a consequence it violates the sum rules which reflect energy-momentum 
and baryon number conservation. 

To restore the partonic picture we use "old fashioned" perturbation theory. That is we 
calculate the g —> cc splitting function P cg in the infinite momentum frame with all three 
partons on-mass-shell. The parton four momenta are shown in Fig. 3. If the momentum of the 
gluon is large, p g 3> kx and m c , then the quark momentum is given by 



zp g + 



ml + k\ 
2zp g 



k T , ZPg 



(16) 



and similarly for k' with z 
splitting in the form 



1 — z and k? — > —k?- We may write the probability of the g — > cc 



dw. 



<■<) 



9 d 2 kT(lk\\ 



(2zp g ) 2 2(1 - z)p g J (AE) 



Sp 



(17) 



contain the normalization factors of the two t 



with colour factor Tr = | and where the [ 
channel and one s channel quark lines shown in Fig. 3. The energy denominators 



AE 



E r 



E n 



ml + k\ 
2z(l - z)p g 



(18) 



5 Of course, to use the parton density for other processes we must calculate the coefficient functions within 
the same framework. 

6 In the Q 2 3> m 2 . domain the evolution in Q 2 is exactly the same for both definitions, however the initial 
conditions of such an evolution, if it were to originate in the threshold region, would be different; compare our 



approach with that of ref. 13 
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play the role of the quark propagators and the numerator 

Q 1 ri rr ( f + m c -ft + m c \ 

= (ml + 4) Z \ + Z ^ z f + < (19) 

where \5^ b is the average over the two transverse polarizations of the (on mass shell) gluon and 
\($ + trie) is the quark density matrix. The factor of 8 in fllTD arises from the sum over two 
polarizations of both the c and c and allows for the t channel parton to be either c or c. 



To identify the splitting function we must rewrite ( |17D in the form 

as dz dQ 2 
2^~~Q 2 



= ^-^P^^Q 2 ) (20) 



where dk\\ = p g dz. The outstanding problem is therefore to determine the scale Q 2 appropriate 
for fcf-. The scale Q 2 should be chosen so that it correctly reproduces the timescale of the 
fluctuations of the gluon into the cc pair, that is 

At ~ AE = # (21 > 
where AE is given by (0). It follows that the appropriate scale would be 

Q 2 = 2 Pg AE = ^|±|. (22) 

Indeed this would be the physically natural scale to adopt. It represents the value of Q 2 
( > Am 2 ) for which the resolution is sufficient to observe the individual c and c partons within 
the short g ^ cc fluctuation time At. 



Unfortunately (|22| ) is not the conventional scale used in the MS factorization scheme in the 
massless quark limit, or rather when Q 2 3> m 2 . To obtain parton densities which correspond 
to the MS scheme for Q 2 ^> m 2 we must take the evolution scale to be[| 

Q 2 = m \ + k 2 T . (23) 

Of course we could use (|22f) as the evolution scale but then we would have to change the 
NLO splitting and coefficient functions. Since the NLO coefficient functions, not only for deep 
inelastic scattering but also for other processes, have been calculated in the MS scheme, it is 
clearly desirable to remain in this scheme. We therefore adopt (|23| ) as the evolution scale. Then 

7 It is connected with the fact that in dimensional regularization in the MS scheme a factor (k^/f^Y is 
introduced into the dk^/k^ integration, see for example ref. Ilffl. Just as the MS and MS scheme coefficient 
functions differ simply by a constant factor, the coefficient functions in the MS scheme and the 'natural' scheme 
based on (E3) differ by a factor which is a function of z. 
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we can use the MS NLO splitting functions P^ for the massless quarks and gluons and for the 
massive quarks for Q 2 m 2 . Moreover for Q 2 m 2 all the NLO coefficient functions have 
the MS form, except for one contribution which we discuss in section 3.2. 



Now using the evolution scale (|23"D we have 



p(o) 

eg 



2T fi 



y + (i 



+ 



2m? 

~Q 2 



-z(l 



m. 



(24) 



Recall that here P cg stands for the splitting into both c(c) and c(c). An analogous result for QED 
may be found in ref. (T^j. The 9 function represents the threshold {k\ = 0) for generating in the 
evolution c and c parton densities, which smoothly tend to the conventional MS distributions 
at high Q 2 . We see that even if at small x we have more than enough energy W to create a cc 
pair, W 2 ~ Q 2 /x ^> Am 2 , then it is possible that the value of Q 2 will be insufficient to resolve 
the c and c pair within the short g <-> cc fluctuation time At, that is when Q 2 < 4m 2 . 

The complete effect of the quark mass in the NLO splitting functions which involve the 
charm quark is contained in fl2"4]). It leads to the following correction to P gg 
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m„ 



T R 5(l-z)^6(Q 2 -m 



(25) 



see (|T5|). 



3. Coefficient functions for deep inelastic charm production 

Just as for light quarks, the contribution of charm to the deep inelastic structure function 
F2 is obtained from a convolution of the parton distributions and the coefficient functions. We 
have 



F 2 C (*,Q 2 ) = ~ £dz- z [c g=c (z,Q 2 ,/i 2 )c(^ 2 ) + C 9 (z,Q 2 ^ 2 )g(- z ,^ 



(26) 



where, due to the quark mass, the coefficient functions have an explicit dependence on Q 2 . The 
charm quark coefficient function in (EB1) has the form 



while for the gluon we have 



C c = Cf) + + (27) 

C,= ^C« + .... (28) 

At NLO accuracy, to which we are working, we need only the coefficient functions that are 
shown explicitly in (|27D and (f2~8|). 

We see that at low scales below partonic threshold, Q 2 < m 2 , where c(x, Q 2 ) = 0, the 
structure function F| is described entirely by jg fusion, that is by the C g (g> g convolution. 
However, we will find that as Q 2 increases from the charm threshold the contribution from the 



S 



7C interaction, C c <8> c, increases rapidly and soon becomes dominant. Of course, as we have 
already mentioned in the introduction, when the number of active flavours increases from 3 to 
4 (as we pass through the threshold region) we must take care to avoid double counting. For 
example, if we were to take the limit in which charm is regarded as a heavy quark, and never 
a parton, then the entire contribution to F2 is 



f = s c r" (8 »- (29) 



We call this fixed (three) flavour approach the photon-gluon fusion (PGF) approximation. From 
the above discussion it might appear that the PGF approximation, which clearly gives the cor- 
rect NLO answer for Q 2 < m 2 ., will dramatically undershoot the true prediction as Q 2 increases 
above the charm resolution threshold. This is not so, since part of the Feynman diagram which 
is responsible for the important C c <8> c parton evolution contribution is contained in C PGF (g) g 
in the PGF approximation (bl. Thus to avoid double counting we will have to subtract this 



contribution from (7 PGF Cg>(7. The consistent treatment of charm mass effects will therefore allow 

9 u 

us to quantify the accuracy of the PGF approximation to F% as a function Q 2 . 

3.1. The charm quark coefficient function for F 2 C 

We must specify the coefficient functions for F% that we introduced in (p6[)-(p8). First the 
LO charm quark coefficient is given by 

C?\z } Q 2 ) =zS(z- 1/(1 + m 2 JQ 2 )) (l + ^pj (30) 

where here z is defined with respect to the charm quark 

( i+ sr <3i) 

The last equality follows directly from the mass-shell condition (x'p + q) 2 = m 2 where x' is the 
fraction of the momentum of the proton that is carried by the struck charm quark, see Fig. 4. 
The final factor in ( j30|) allows for the Fl component of F2 = Ft + Ft where 

a L /a T = Am 2 JQ 2 . (32) 



Inserting Cf ] of © into © gives a contribution to F| (x, Q 2 ) proportional to xc(x', Q 



2\ 



where here the true scale is fx 2 = Q 2 . In fact at NLO all the m c 7^ effects in the charm quark 
coefficient function occur in Cjp\ Indeed we justify in section 3.3 that at NLO we may simply 
use the massless quark expression for the coefficient . 

3.2. The gluon coefficient function for F^ 

We may write the gluon coefficient function, defined in (p8[), in the form 



Cf = CJ GF - AC g (33) 



where the PGF expression for F 2 is [p"9 



+ 



z 2 + (1-z) 2 



k(l - z) 



z 1 — 3z 



5WI 



c _2 



Q' 1 



hi 



1 + 

1-/3 



0> #(Q 2 (--1) - 4m^).(34) 



/3 is the velocity of one of the charm quarks in the photon- gluon centre-of-mass frame 

4m 2 . z 



3 1 



1 - 



Q 2 {l-z) 



(35) 



The 9 function in (0), 9{W 2 — 4m 2 ), represents the cc production threshold, where W is the 
cm. energy. Its presence guarantees (3 2 > 0. The AC S term in (|33| ) is necessary to avoid the 
double counting of the graph that we have already used to compute Pj$, see section 2.3. That 
is we must subtract from C FGF the term P^ Cg> CjP^ that we already include in the parton 
evolution up to Q 2 . The z variable in the gluon coefficient functions is defined with respect to 
the gluon momentum fraction x g , 

(36) 



X 



.1 n 



Z Z 



where z' = x'/x g , see Fig. 4. Thus the explicit form of the subtraction term is 

dzo 
z 



&C g (z,Q 2 ) = { ^CjP\z ,Q 2 ) f d\nQ' 2 P cg (z',Q' 2 ) 

J Zn JQi 



1 + 



4m 2 c 



( Q d\xiQ' 2 P cg (z',Q' 2 ). 
Jo 2 . 

^ mm 



(37) 



The lower limit of integration is given by the "partonic" 9 function which is hidden in P cg (z', Q' 2 ^ 



that is 



Q 



2 

min 



m 2 . 



The integration in ([J7]) may be readily performed to give 



AC g (z,Q 2 



[z 12 + (1 - z') 2 } In 



m 



4mV 
+ ^ 



(38) 



(39) 



where we require Q > Q 2 nin and where z' = z/z = (1 + m 2 /Q )z 



It is interesting to consider the m 2 — > limits of C g GF and AC g . We have 



^PGF 
°9 



{z 2 + (1-z) 2 } lnf— ^) + 8z(l 

\ Z m c ) 



Z)-1 



(40) 



as m c — > 0, which differs from the exact m c = coefficient by the presence of Q 2 /m 2 c in 



the argument of the logarithm. However, from (|39"D we see that 

AC g - {z 2 + (l-z) 2 } ln(^ +2,(1-,) 



(41) 
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as m c — > 0, which removes the \n(Q 2 /m 2 ) term in = Cj GF — AC r 

The difference C PGF — AC g is exactly as in the conventional (massless) MS scheme apart 
from the final term 2z(l — z) in AC g . This discrepancy is due to the different order in which the 
limits are taken in the calculation of the coefficient function. In the massless case we first take 
m q — > and then e — * 0, which is equivalent to taking the infrared cut-off A > m q , whereas here 
we have first to take e — > and then consider the Q 2 3> m 2 limit. That is we can only neglect^ 
m 2 in P(g> of (|24D if Q 2 ^> m 2 ,, whereas we need to evaluate P cg in fl37|) down to Q 2 = Q^in = m c- 



3.3. Choice of scale 

We now come to the choice of the scale /x 2 in fl26|). First we consider the convolution involving 
the gluon. The only dependence on the scale /i 2 in the coefficient function C g is in the argument 
of as in Cj} 1 ' has no dependence on /i 2 since all the collinear singularities are regularized 
by the heavy quark mass provided that /i 2 < m 2 . Nevertheless we have to choose the scale /z 2 for 
as and the gluon distribution. Variation of the scale induces only NNLO contributions. There 
is as yet no complete calculation of the NNLO contributions^ (in our framework) applicable for 
all Q 2 which would introduce terms compensating the variation with scale. We must therefore 
attempt to identify the 'natural' scale for the process. We have already mentioned that the 
natural scale for the charm convolution is fi 2 = Q 2 . 

Due to the different way that the scales enter, the as(fi 2 )AC g <8> g(x/z,fi 2 ) term does not 
exactly subtract the LO charm contribution which comes from the convolution CfP> <S> c(x, Q 2 ). 
At first loop level the latter term is of the form 

c ' o) Cv asiq2)p ^ ){z)g & q2 ) • (42) 

Here we also take the natural choice]^] of scale fi 2 = Q 2 in the AC g term. We see that over 
the whole range of integration in (|i^) we then have q 2 < fi 2 and as(q 2 ) > as(fi 2 ). In other 
words the subtracted value of the AC g term is a little less than needed, leading to a lack of 
smoothness in F 2 C near threshold. To diminish this effect we could reduce the scale fi 2 in the 
AC g term by taking, say, [i 2 = 5Q 2 with 5 < 1. From the formal point of view the choice of 
fj, 2 should not matter. In a NLO analysis it only induces changes at NNLO. However m 2 is not 
so large, and some of the deep inelastic data used in our fit will be at sufficiently low values of 
Q 2 that we sample scales /j, 2 > m 2 . In this Q 2 domain the analysis does have sensitivity to the 
choice of scale, showing the need for the NNLO formulation. 

Now let us return to the charm quark coefficient function and, in particular, explain why the 
massless approximation is sufficient for at NLO. After the subtraction of the LO contribu- 
tion (in an analogous way to (|33| ) for the gluon coefficient function) the remaining coefficient 

8 If we were to neglect rr? c in ( pi ) for all Q 2 then the additional 2z{\ — z) contribution in AC g would disappear. 
9 However for Q 2 » m 2 see |l3| . 

10 Of course we do not want [i 2 < m 2 and so, in the analysis described in section 4, we take [i 2 = max{Q 2 , m 2 }. 
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contains no InQ 2 terms. Recall that in the presence of a heavy quark mass we lose a 
logarithm (see ( Jl2| ) and (0)). Therefore at NLO we need only consider m c 7^ in the LO 
diagrams. Similarly at NNLO we only require m c 7^ in the NLO graphs and so on. Therefore 
the (non-logarithmic 0(as)) contribution is only needed when Q 2 ^> m 2 , as was discussed 
in section 2.2 after Eq. fll3l) . For Q 2 < m 2 the charm density c(x, Q 2 ) = 0, while for Q 2 > m 2 
two powers of InQ 2 are lost (one in P cg and one in C^) and so this region contributes only at 
NNLO. Thus we can set m c = in C^' . Of course we must use the same definition of the scale 
Q 2 , (|23|) , as for massless evolution in the MS scheme. 



3.4. The resolution of charm 



Although we now have a definite framework which enables us to incorporate a m c 7^ charm 
parton into a parton analysis, we immediately encounter a problem when we confront the data. 



The charm density rises rapidly as we evolve up from the threshold Q 2 



mt 



(required by the 



MS scale (|23|)) such that in the region Q 2 > m 2 it is in conflict with the data. The reason is 
clear. The physically reasonable scale at which the charm parton may be resolved is given by 



(2~2~) with a threshold at Q 2 = Am 2 . In order to implement this behaviour we introduce a factor 



./ 



Am 2 

"o 1 



1- 



Am 2 

"o 1 



(43) 



into the charm coefficient function C c of (|2T| ) and (|30| ) which, via (|37|), then feeds through into 



AC g . 



At first sight the introduction of the 'ad hoc' factor of / appears to modify even the LO 
result, which in symbolic form now may be written 



F 2 c (LO) = fCW <g> c . 



(44) 



This is not so. In fact we will see that the modifications due to the introduction of / only enter 
at NNLO. At LO we have strong ordering in transverse momenta. LO contributions therefore 
only occur for Q 2 3> m 2 , where / — > 1. This reflects the fact that mass effects correspond to the 
loss of a factor of InQ 2 and only contribute at the next order of as- So the NLO contribution 
is changed, but only by an amount ~ as ln(4m 2 /m 2 ) ~ const, as coming from evolution over 
the limited interval m 2 < Q 2 ~ Am 2 . However even this contribution is cancelled by the AC g 
coefficient function. To see this we inspect the NLO form 

F 2 c (LO + NLO) = / (Cf ) + a s C^) ® c + a s C^ ® g . (45) 

To obtain the gluon coefficient function recall that we evaluated jg — * cc at 0(as), which 
we denoted by asCg GF , and then subtracted the LO part (~ a^lnQ 2 ), which was already 
included, 



F 2 c (LO) = fC®®c 

= f Cf ® a s In Q 2 P^ ®g = a s AC^ ® g 



(46) 
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where it is sufficient to use the LO expression for the charm density. Thus 

= a s C^ GF -fC^ ®a s \nQ 2 P^. (47) 

The last term of ( f47"D exactly cancels the first term of (flop , so 

F 2 c (LO + NLO) = a s C PGF ® g + fagC® <g> c . (48) 

Thus the introduction of the ad hoc factor / gives rise to a modification (1 — f)asC^ ® c 
which only enters at NNLO. Indeed, in the appendix we show precisely how the modification 
is cancelled when working to NNLO. 

Thus, in summary, in the NLO global parton analysis of section 4 we make the replacement 

C< 0) -> fC® (49) 
in (j30D , and similarly for Cj^\ Hence (|39"D becomes 

AC g -> fAC 9 , (50) 



where / is given by fl4"3"|). After the introduction of the factor /, only the PGF contribution 
survives in the region Q 2 < Am? c below the resolution threshold, even though we have a non-zero 
charm parton density for Q 2 > m 2 . As we evolve above the resolution threshold Q 2 = 4m 2 the 
charm parton component of F 2 rapidly becomes important. 



4. Charm as a parton in a global analysis 

The measurements of F 2 at HERA have become much more precise with errors as small as 
±3% or less. Moreover, since the charm component _F 2 C of F 2 is about 0.25 in the HERA regime 
it is important to improve the treatment of charm in the analysis of deep inelastic scattering 
data. This was the objective of sections 2 and 3 above. The new formalism incorporates 
the heavy quark masses in the parton evolution equations and allows a determination of the 
(universal) charm and bottom quark densities. Indeed we can predict c(x, Q 2 ) and b(x, Q 2 ), as 
well as the charm and bottom components of F 2 , directly from a knowledge of the gluon and 
other quark densities. There are no free parameters, although the results do depend on the 
values of m c and mj, and, as usual, on the truncation of the perturbation expansion. As in 
previous analyses, we work to NLO. 

The new framework is a significant advance on the existing treatment of charm in deep 
inelastic scattering. Recall that two different types of approach are used at present. In the 
first, charm is set to zero below some scale (c(x, Q 2 ) = for Q 2 < fi 2 ) and for Q 2 > /i 2 the charm 
distribution is evolved assuming that m c = 0. Although this procedure is clearly inaccurate 
in the cc threshold region, the parameter \x is chosen so that the fixed-target F 2 data are well 
described. Secondly, we have the PGF approach |10|, [11] based on the calculation of 7*g — > cc 



13 



with the correct kinematics, but in which c is not treated as a parton. As we have seen, this 
gives the correct description of F 2 for Q 2 < m 2 c and should remain a reasonable approximation 
to F 2 for Q 2 > m 2 . However, the PGF model will inevitably break down at larger Q 2 values 
when charm can no longer be treated as a non-partonic heavy object and when it begins to 
evolve more like the lighter components of the quark sea. 

Actually there exists in the literature a third, hybrid, approach [|14]]. Charm is treated as 
a new massless parton above Q 2 = m\. That is the ml effects are neglected in the splitting 
functions, although they are included in the coefficient functions to NLO. This is not quite 
correct since the neglected m 2 effects would give NLO contributions during the evolution. 

Before we present our predictions for c(x, Q 2 ) and b(x, Q 2 ), we perform a NLO global analysis 
of deep inelastic and related data which incorporates the m q ^ parton evolution procedure 
that we introduced in sections 2 and 3. This may be regarded as a small refinement of the 
global analysis determination of the gluon and light quark densities of ref . || , but it does allow 
the gluon (and other parton) distributions to readjust themselves to accommodate the new 
treatment of c(x, Q 2 ). Recall that the heavy quark distributions, c(x,Q 2 ) and b(x,Q 2 ), do not 
contain any free parameters apart, of course, from m c and m&. Motivated by QCD sum rules, 
we take m c = 1.35 GeV and m& = 4.3 GeV [M. We show the effects of varying the value of m c 
when we discuss the description of F 2 . In fact we find that the overall description of the data 
(and in particular of F2 in the HERA regime) improves compared to our previous analyses ||. 
The only change to the data set that we use is the addition of the final NMC data [2]J for F 2 . 

We shall present full details of the new global analysis []] in a future paper in which we 
will discuss the improvements of the deep inelastic data and their implications. However, in 
Table 1 we illustrate the quality of the new fit relative to our previous fit that best described the 
HERA data, MRS(R2) |§. We see that despite now having a prescribed charm distribution, the 
quality of the new fit is comparable to or actually slightly better (particularly for the small-x 
F 2 measurements) than that of the previous analysis. 

The HERA data lie in the region where F 2 /F 2 is largest and there is clear improvement 
in the new fit for these data. The value of as resulting from the new fit is as(M§) = 0.118, 
intermediate to the values 0.113 and 0.120 of MRS(Rl) and (R2) and the lower \ 2 f° r the 
BCDMS data in the Table is due to this. Our prescription for cts(Q 2 ) across charm and 
bottom thresholds is to match the values at Q 2 = m 2 , and again at Q 2 = m 2 . Thus we define 

a5(4)(Q 2 ) = "s(Q 2 ,4) (51) 

and take, for 5 flavours, 

« 5 (5)W 2 ) = as\Q 2 , 5) + as\ml 4) - a s \m 2 b , 5), (52) 



11 The FORTRAN code for this set of partons, MRRS, together with the code for computing each flavour 
component to F\, F2 and Fl is available by electronic mail from W.J.Stirling@durham.ac.uk, or directly from 



frttp : //durpdg . dur . ac . uk/HEPDATA/MR£ 
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Experiment 


7^ data 


MRRS 


x 2 

A. 

MRS(R 2 ) 


HI 


193 


133 


149 


ZEUS F 2 ep 


204 


290 


308 


BCDMS 


174 


271 


320 


NMC F$ p 


130 


145 


134 


NMC F 2 M 


130 


119 


98 


E665 F 2 W 


53 


60 


62 


E665 Ff 


53 


54 


60 


SLAC F 2 ep 


70 


96 


95 



Table 1: \ 2 values for some of the data |§ g, |T|, [|||23], |_ 
larger \ 2 values for the E665 points |24j than those quoted in ref. 



used in the global fit. Note the 
- these result from our 



correcting our previous incorrect treatment of the E665 experimental errors. 



while for 3 flavours we have 

«5(3) (Q 2 ) = ®s\Q 2 , 3) + ag\m 2 c , 4) - c^Vc, 3). (53) 

In Fig. 5 we show the flavour decomposition of the S6cl clS db function of Q 2 for two different 
values of x. Recall that there are now no input parameters for the heavy quark distributions, 
c(x,Q 2 ) and b(x,Q 2 ), and that they are determined in terms of the gluon (and other parton) 
distributions. 

We show the description of both the fixed target and HERA data for F 2 in the next section. 
The charm data are not used in the global fit. However, when they become more precise these 
data should be included as they will provide a significant extra constraint on the gluon distri- 
bution. The gluon density from the new fit compares very closely with that of MRS(R2). The 
new gluon is more 'valence-like' at Ql = 1 GeV 2 , but for Q 2 > 2 GeV 2 both gluon distributions 
rise at small x and become increasingly similar as Q 2 continues to increase. 



5. The structure of F 2 C 

Fig. 6 shows the partonic decomposition of F 2 C as given by (p6|), which may be written in 
the symbolic form 

F 2 C = C c ® c + C g ® g. (54) 

The gluonic component gives the total production below the charm resolution threshold, Q 2 < 
4m 2 . However, the component driven by the charm distribution rises rapidly above threshold 
and becomes dominant at larger Q 2 . We also show for comparison the photon-gluon fusion 
prediction C PGF ® g. The PGF model and our prediction are identical below threshold, Q 2 < 
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4m 2 . Above threshold we see that the rapid onset of the charm parton component C c ® c 
is largely balanced by the subtraction AC 9 from the PGF result. Let us discuss in turn the 
behaviour of _F 2 C near the charm threshold and then at large Q 2 . 

The lack of smoothness of F% apparent in Fig. 6 in the charm threshold region is due to the 
mismatch of the subtraction term asAC g (g) g with C^> <S> c = ® ocsPj^' <8> g with different 
scales of as and g in the two terms, see section 3.3. In the AC g term the scale is fx 2 (where 
we have taken the natural choice /x 2 = Q 2 , or rather /i 2 = m&x{Q 2 ,m 2 }), whereas as and g 
in the second term are evaluated at scales varying over the convolution interval m 2 to Q 2 . Of 
course we could have reduced the mismatch by choosing a smaller scale /z 2 , more representative 
of the m 2 to Q 2 integration interval. But formally in a NLO analysis the choice of scale /i 2 
(and the mismatch) should not matter. It only gives contributions at NNLO. However in the 
charm threshold region as(fi 2 ) is relatively large and we are sensitive to the choice of fi 2 . If 
the NNLO formalism were available the behaviour of F£(x, Q 2 ) in the charm threshold region 
would be more stable under variations of /i 2 , and would have a smoother form in Q 2 . 

As expected these problems evaporate at larger values of Q 2 . Away from the charm threshold 
region (Q 2 > 20 GeV 2 ) the predictions for F 2 C for different fi 2 rapidly approach each other as 
Q 2 increases and become insensitive to the choice of scale fi 2 . The effects of the evolution of the 
charm density are evident. A measure of the effect is the difference between the prediction of 
F 2 (continuous curves in Fig. 6) and that obtained in the PGF model (dot-dashed curves). By 
Q 2 = 100 GeV 2 , for example, for x = 0.05 (0.005) the improved description, in which charm is 
treated as a parton, lies some 75% (30%) above the PGF model. 

The comparisons of the predictions for F 2 with the EMC and the HERA data are shown in 
Fig. 7. The overall agreement over quite an extensive range of x and Q 2 is good. The dotted 
and dashed curves in Fig. 7 show the effect of taking m c = 1.2 and 1.5 GeV respectively, rather 
than the central value, m c = 1.35 GeV, which we use throughout this paper. 

Fig. 8 shows the fraction of charm deep inelastic events as a function of Q 2 for selected 
values of x. The strong production of charm at HERA is evident; moreover we see a sensitive 
dependence on x and Q 2 . If a significant fraction of the numerous charm events can be cleanly 
isolated in the experiments at HERA then the resulting precision measurement of F 2 , coupled 
with the measurement of F 2 , will provide a powerful double constraint on the gluon distribu- 
tion, as well as offering a stringent scheme independent test of QCD along the lines of that 



using F 2 and F L proposed by Catani [25 



6. Predictions for F[ 

We may also use the new formalism which incorporates the quark mass to calculate the 
charm component F[ of the longitudinal structure function. We use expressions that are iden- 
tical to fl2HD-(]2"8D and (|3"3] ) but with the coefficient functions C q=c and C g that are appropriate 
to F£. For the quark coefficient we have 

C7?> = ^z5(z- (l+m 2 /^)- 1 ), (55) 
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whereas for C^ 1 ' we may use the massless quark expression, since we are working to NLO 
accuracy. For the gluon coefficient for F[ we have 



U 9 



where 



with Q 2 
term is 



C. 



PGF, 



z,Q 2 



Af3 z(l 



-^r In 



1+/3 



Q 2 l 

> Am 2 c z/(1 — z), where the quark velocity (3 is given by 







(56) 



(57) 



Here the subtraction 



&C g (z,Q 2 



Am 2 



(58) 



where [. . . ] is the expression in the square brackets in ([39D . For AC g to be non-zero we require 
Q 2 > Qmin' where Q^in i s given by (^). Just as for the coefficient functions for F 2 , we also 
include the factor / of (|3|) in and AC g . 

In Fig. 9 we present the predictions for Fl in terms of the ratio R c = F[/F^. Due to the 
factor 4m 2 ./ Q 2 in the coefficient function of the LO charm component given in (|5l|), we have 
a pronounced peak in R c just above the resolution threshold, Q 2 = Am 2 . In this region R c 
is sensitive to the precise choice of the scale /i 2 . As expected R c decreases as Q 2 increases, 
as well as becoming more stable to changes of scale. The NLO gluonic component gives a 
smaller value of R c than the charm component. Hence the peak is more pronounced at larger x 
when the gluonic component is less important. We also show in Fig. 9 the values of R = Fl/Ft- 



7. Conclusions 

We have determined the charm and bottom quark densities of the proton taking into ac- 
count the effects of their non-zero mass. In particular we have presented a formalism which 
incorporates m c and m^ into the Altarelli-Parisi splitting functions and in the coefficient func- 
tions in a consistent way. We can therefore evolve up in Q 2 taking proper account of the heavy 
quark thresholds. At NLO accuracy we show that the main effect of the quark mass is in the 
splitting function Pf® (or Pf^)- 

We showed that the threshold for the charm density, c(x, Q 2 ), occurs at Q 2 = m 2 . On the 
other hand we know that the threshold for deep inelastic cc production is given by W 2 = Am 2 , 
or equivalently Q 2 = Am 2 (l — x)/x, which for small x occurs below the partonic threshold 
Q 2 = m 2 . This apparent contradiction has a simple explanation. In the region Q 2 < Am 2 we 
find that Q 2 is too small to allow sufficient time to observe the g — > cc fluctuations which occur 
within the proton. Here the photon-gluon fusion mechanism, j*9 ~^ c c, gives the complete 
answer. For evolution above the partonic resolution threshold the structure of F 2 C is more 
interesting. The charm component 7*c — > c with a spectator c quark (or vice-versa with c <-> c) 
increases rapidly and soon exceeds the gluonic contribution 7*g — > cc which only enters at NLO. 
In the partonic description the LO part of the gluon now has the structure (g — > cc) <8> (7*c — > c). 
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To avoid double counting we must therefore subtract this LO contribution of the gluon and 
keep only the part coming from Cjj 1 '. 

In addition to its importance in determining the charm quark density c(x,Q 2 ), the correct 
formulation of charm mass effects in evolution has become essential in order to obtain an 
accurate description of F2 in the HERA domain. The reasons are that the charm component of 
F2 is appreciable (F2/F2 ~ 0.25 for x ~ 0.001 and Q 2 ~ 25 GeV 2 ) and that the measurements 
of F 2 at HERA are now much more precise. 

In summary, in this paper we have shown how to treat charm as a parton for all values of 
Q 2 . The new NLO partonic formulation, which incorporates m c 7^ effects, has the following 
important features. 

(i) The charm distribution contains no free parameters, except m c . 

(ii) The partons are universal (that is they can be used in the NLO description of all hard 
scattering processes initiated by protons). 

(hi) the splitting and coefficient functions coincide with those of the (massless) MS scheme 
for Q 2 3> m 2 (with the one exception discussed in section 3.2). 

(iv) The momentum and flavour sum rules are conserved. 

(v) There is a definite prescription to enable the formulation to be extended to include m c 7^ 
effects at NNLO and higher orders. 

(vi) The new framework, in which the charm density is defined in terms of a leading lnQ 2 
decomposition of the Feynman diagrams retaining the full mass effects, is applicable in 
the important threshold^ region Q 2 > m 2 . 

Finally, we have used the new prescription to perform a global analysis of deep inelastic and 
related hard scattering data and generated charm and bottom quark densities. The analysis 
predicts the values of F 2 C (and -F|). We find that the predictions for F^ show some sensitivity to 
NNLO effects in the charm threshold region (Q 2 ~ m 2 ), but become increasingly stable as Q 2 
increases above about 20 GeV 2 . We find good agreement with the EMC and HI measurements 
of F 2 C . These data, which span a wide range of (x,Q 2 ), were not used in the global analysis. 
Clearly as the experimental precision increases they should be included, and will impose a 
valuable additional constraint in the determination of the parton densities, and of the gluon in 
particular. 
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12 The latter region is inaccessible to the RG approach to the charm density. Indeed it is not clear how 
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formulation reduces to the conventional RG massless parton approach. 
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Appendix 

Here we demonstrate how, if we were to work at NNLO, the factor / of (fS|) contributes 
only at NNNLO. If, for simplicity, we neglect the light quarks, then we have in analogy to 



F 2 % . . + NNLO) = / (C<°> + asC^ + a|Cf } ) ® c + (a s C^ + o? s Cf>) ® g . (Al) 

C^p is given by the 0(a s ) expression for the photon-gluon cross section, a|C^ 2 ) PGF minus 
the a|ln 2 Q 2 and a 2 s \nQ 2 contributions which are already generated within the LO+NLO 
formalism. That is 

a 2 C^ = afCf PGF - ACf (A2) 



with 



eg ^ 1 99 ^ 1 cc ^ 1 eg J ^ y ^ "5 111 V 1 eg 

+fa s CP ® a s In Q 2 Pj°) ®g + a s C^ ® a s In Q 2 P^ ® g . (A3) 

Inserting (A2) into (Al) and cancelling terms, we find that the residual 0{a 2 s ) part of _F 2 C is 

F 2 c (NNLO) = a 2 s C^ PGF ® g + fa 2 s C^ ® c , (A4) 

in analogy to ([ISP . The modification (1 — f)a 2 s C^ ® c due to the introduction of the ad hoc 
factor / is now of NNNLO. 
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Figure Captions 



Fig. 1 Part of the parton chain occurring in the description of deep inelastic scattering which 
contains the g — > cc transition. 

Fig. 2 An example of a "block" diagram along the parton chain, which gives NNLO charm mass 
effects if the two s channel charm quarks have comparable transverse momenta. Then 
the charm mass should be retained for all the quark lines that are shown. 

Fig. 3 The diagram used to calculate the charm mass effects in P} g '- 

Fig. 4 The variables used in the discussion of the coefficient functions C q=c (z, Q 2 ) and C g (z, Q 2 ). 
For the charm quark function the variable z = x/x', while for the gluon function z = x/x g , 
see eq. ( pip and ( ]3T)| ) respectively; x is the usual Bjorken x = Q 2 /2p.q. 

Fig. 5 The flavour decomposition of the quark sea distribution of the proton as a function of Q 2 
at two values of x. The total sea is given by S = 2(u + d + s + c + b) . 

Fig. 6 The partonic decomposition of F 2 C as a function of Q 2 for x = 0.05 and x = 0.005. For 
Q 2 < 4m 2 there is only the contribution from C g = C FGF . For larger Q 2 , C g = C g — AC g 
and the total F£ is the sum of this contribution and that from C c . 

Fig. 7 The description of the EMC and HERA measurements of F%. The solid line corresponds 
to our new fit with m c = 1.35 GeV. The dashed and dotted lines correspond to taking 
m c = 1.5 and 1.2 GeV respectively, with all other parameters unchanged. 

Fig. 8 The ratios F2/F2 and F2/F2 at fixed values of Q 2 resulting from the new global fit (in 
which we take m c = 1.35 GeV and rrib = 4.3 GeV). The experimental data point shows 
the estimate from ref. for F^jF-i in the HERA kinematic region, 10 GeV 2 < Q 2 < 
100 GeV 2 . 

Fig. 9 The predictions for R c = F[/F£ and R = Fl/Ft as a function of Q 2 for x = 0.0005 and 
x = 0.05. 
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